Boson Hubbard model with weakly coupled Fermions 
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Using an imaginary-time path integral approach, we develop the perturbation theory suited to the 
boson Hubbard model, and apply it to calculate the effects of a dilute gas of spin-polarized fermions 
weakly interacting with the bosons. The full theory captures both the static and the dynamic effects 
of the fermions on the generic superfluid-insulator phase diagram. We find that, in a homogenous 
system described by a single-band boson Hubbard Hamiltonian, the intrinsic perturbative effect of 
the fermions is to generically suppress the insulating lobes and to enhance the superfluid phase. 
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Introduction. The boson Hubbard model has long pro- 
vided the paradigm for studying one of the simplest quan- 
tum phase transitions (QPT), the superfluid to insulator 
transition (SIT) in a dilute gas of bosons. Most satis- 
factorily, recent experiments 1,2 using ultra cold bosonic 
atoms confined to an optical lattice, which mimics the 
boson Hubbard model in a custom setting, demonstrated 
the existence of the SIT in a pristine, disorder free, boson- 
only system. By varying the effective jj of the ultra cold 
atoms in optical lattices, where t is the boson nearest- 
neighbor hopping parameter and U is the on-site boson- 
boson repulsion, the researchers demonstrated the exis- 
tence of the Mott-insulating (small jj) and the superfluid 
(large jj) states in the time-of- flight experiments ~£ At 
some intervening value of jj , then, there should be a QPT 
separating the two states^ 

An important theoretical question, which has received 
wide attention££i£&i2*ii in light of the recent experi- 
ments in the Bose-Fermi mixtures , 12 i 13 is what happens 
to the insulating and the superfluid phases when fermions 
are introduced to the bare boson Hubbard model. In the 
case of the bosons weakly interacting with spin-polarized 
fermions, which are away from half- filling, this question 
can be addressed analytically. While some of the ear- 
lier studies^^ concluded that the region occupied by the 
superfluid phase in the phase diagram is enhanced by 
fermions, more recent onesii concluded that the oppo- 
site is true because of an effect akin to the fermionic 
orthogonality catastrophe due to the dynamic effects. In 
this Communication, we address this question by devel- 
oping a rigorous perturbation theory suited to the single- 
band boson Hubbard model, which captures both the 
static and the dynamic effects mediated by the fermions. 
Our conclusion is that, in a homogenous, single-band sys- 
tem and in the absence of loss of cooling due to adding 
fermions, the fermions intrinsically shrink the area occu- 
pied by the Mott insulating lobes (Fig. [T]), thus generi- 
cally enhancing the superfluid region. The overall effect 
is qualitatively in the same direction as in the effects 
of Ohmic dissipation in enhancing the superconducting 



phase coherence in Josephson junction arrays^ or in 
granular superconductors^ Even though in the current 
experiment a 12 ! 13 the Bose-Fermi interaction strength is 
not in the perturbative regime, it is possible to tune this 
coupling and bring it to the perturbative regime— Thus, 
our predictions can be tested experimentally. Further- 
more, in light of our present analytical results (and the re- 
sults of Refs. [ETflol]) it seems likely that the observed loss 
of superfluid coherence by adding fermion o 12 ' 13 should 
be attributed to the external factors, such as heating 17 
and self-trapping of the bosons and fermionsJ^ Hence, 
experiments which can avoid such effects (e.g., shallower 
lattices and lower boson filling factor have reduced boson 
self trapping due to fermions^) are necessary to see the 
intrinsic effect - enhancement of the superfluidity - due 
to the fermions. We stress that the perturbation the- 
ory of the boson Hubbard model we develop, which de- 
viates from the standard machinery^ applicable to the 
free bosons, should have other important applications, 
e.g., the phase diagram of the boson Hubbard model in 
the presence of coupling to a dissipative Ohmic bath 2 ^ or 
a second boson species. In general, our method, specifi- 
cally Eqs. (|I2[ IT3l [T4]) . can be taken over in any problem 
where the Green's function of the boson Hubbard model 
has to be calculated in perturbation theory. 

The model and the results. We consider a mixture of 
bosonic and spin-polarized fermionic atoms in an optical 
lattice. The Hamiltonian of the Bose-Fermi system is 
written as H = Hb + Hp + Hbf, with 



ff^Elf^- 1 )-^)-* E { b l b J+ H - c -)> (!) 

i ^ ' <ij> 

H F = -t F -M^I>K (2) 

<ij> i 

H FB = U FB ni(c\ci - n° Fi ). (3) 



Here cj and b\ are the fermion and the boson creation op- 
erators on site i, hi = b\bi is the boson density operator, 
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FIG. 1: (Color online) a) Main panel: Phase Boundary of 
the boson Hubbard model with and without the fermions 
for the boson density no = 1. Solid line describes the 
insulator-superfluid phase boundary without the fermions. 
The dashed line corresponds to the same phase boundary with 
the fermions present. The dash-dot line denotes the phase 
boundary in the static approximation. The regions near the 
degeneracy points (integer (i/U) are implicitly excluded from 



this figure 22 

Inset: The dependence of the function R(y) on its argument 



Here we used -r^r- 



0.1 and ^ =0.15. b) 



U > 0(Ufb) describes the on-site boson-boson (boson- 
fermion) interaction, t(tp) corresponds to the hopping 
matrix element for the bosons (fermions), n Fi is the aver- 
age density of the fermions, and /i = Ho — UpBripi and 
are chemical potentials for boson and fermions, respec- 
tively. Here no is the boson chemical potential without 
the fermions. 

The partition function of the bare model (without the 
fermions) can be written in terms of an imaginary-time 
path integral over a complex scalar field i/>(x, t)^ where 
r is the imaginary time. The action in terms of ^(x, r) 
takes the form of a 4> 4 theory, see Eq. © . In this descrip- 
tion, the details of the bare Hamiltonian are hidden in the 
coefficients of the various terms of the action. For exam- 
ple, the coefficient, r, of the term |?/>(x, t)| 2 (see below) 
is determined by the Green's function, (T T 6i(r)6](0)), 
of the bosons^, where (...) denotes average with re- 
spect to the on-site part of the boson Hubbard Hamil- 
tonian. In mean field theory, r = gives the locus of 
the insulator (r > 0, (i/j(r, r)) = 0) to the superfluid 
(r < 0, (tp(r, t)) ^ 0) QPT, revealing the Mott insulating 
lobes in the phase diagram.—^ 

With fermions, a similar description of the partition 
function still holds, but now the boson Green's function 
must incorporate the perturbative effects of the boson- 
boson interaction mediated by the fermions. We stress 
that this mediated interaction is manifestly non-local 
in both space and time. Therefore, it is not obvious 
that this problem can be treated in an effective Weiss- 
like single-site theory as done in Ref. The pertur- 
bative corrections to the boson Green's function can- 



not be calculated by using the standard diagrammatic 
machinery^ either, because the bare Hamiltonian is an 
interacting one and the interaction U has to be treated 
non-perturbatively. We solve this problem by noting that 
we can still calculate the needed correlation functions 
exactly by making use of the eigenstates of the number 
operators {rii}. A modified linked-cluster theorem still 
holds which gets rid of all the divergences encountered 
in the perturbation theory. The locus of the equation, 
r' = 0, where r' includes the perturbative corrections to 
the boson Green's function, provides the phase boundary 
between the superfluid and the insulating states. Our 
central result for the phase boundary is shown in Fig.l. 
Below we give a summary of the methods and the calcu- 
lations used to arrive at the results. The details of the 
calculations will be given elsewhere^ 

Summary of the methods. To the lowest order in 
Uf b , the effect of the fermions on the constituent bosons 
is a trivial shift of the boson chemical potential fi 
1 1 ; ; — UFBripi- All the non-trivial effects appear in the 
second order in Ufb- By integrating out the fermions, 
the imaginary-time partition function becomes (we as- 
sume here zero temperature T — > 0) 



Z = J ®b*®bi exp (S eS [b*, bi\) (4) 
!S [b*,h] = £ dr Kdrbi + Hb^J (5) 

-} y I dri dT 2 n i (T 1 )M ij (T 1 -T2)n j (T 2 ). 
,■„■ Jo Jo 
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In the second order in U fb, the integral over the fermion 
degrees of freedom gives rise to an effective non-local 
density-density interaction for the bosons with the func- 
tion Mij{r\ - r 2 ) being, 

U 2 

Mijin - r 2 ) = -p (An F i(n)An Fj (T2)) . (6) 

In the frequency and momentum domain, M q (£l n ) is pro- 
portional to the fermion polarization function, and in 2D 
is given by, 



M q (n n ) 



u 2 

U FB 

2A 



(7) 



Here, v n = £l n /AEp and k = q/2kp, with Ep and hp be- 
ing the Fermi energy and the Fermi momentum, respec- 
tively. A is the fermion mean level-spacing, A = 1/vpV, 
with vf the density of states at the Fermi level and V the 
volume of the unit cell. Equation 1(7} is valid for k < 1. 
Here, for simplicity, we consider a 2D system. However, 
our qualitative conclusions hold for the 3D case as well 21 . 

Using the Hubbard-Stratonovich transformation with 
a complex scalar field ipi (r) , we integrate out the bosonic 
fields to write, Z = Z J Qip^if;* exp(-S[ipi,ip*}), where 
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the action S[ipi,ip*] is given by, 
Jo 



(8) 



ln^exp 
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Here the matrix elements of the symmetric matrix, uiij , 
are equal to t for the nearest neighbors and zero oth- 
erwise. The expectation value in Eq. ((8]) is taken with 
respect to the action S e g [b*, &<] (with t = 0). By expand- 
ing i5[y>, V 1 *] up to the fourth power of the field ip, and 
taking the continuum limit, we arrive at the action of an 
effective complex ip 4 field theory, 



S[W]= 



dxlci^ii — +c 2 



9^ 



-c|V^| 




with x = {r, r}. The coupling constants Ci, C2, c, r, u are 
given by the correlation functions of the boson Hubbard 
model with t = 0. In mean field theory, the phase bound- 
ary between the superfluid and insulating states can be 
obtained by setting the coefficient r to zero: 



1 & 

- + / drGi r 



= 0, 



(10) 



where G*(r) = -(T r 6 i (r)6 l t (0)) is the single-site boson 
Green's function, which, in the presence of the fermions, 
should include the effective fermion-mediated density- 
density interaction. Without the fermions, this Green's 
function is given by,* 



Gi(iw n )-- 



(n + l) 



n 



-8E h 



(11) 



where 8E P and 8Eh are particle and hole excitation ener- 
gies: SE P = UriQ — /i and 8Eh = ±i — U(tiq — 1), and no 
is the number of bosons per site minimizing the ground 
state energy. Thus, the problem is now reduced to the 
calculation of the on-site full boson Green's function by 
computing the corrections to Eq. (fTTj) . As we show be- 
low, this can be done perturbatively in Ufb- 

The calculation of the perturbative corrections to the 
boson Green's function is non-trivial because the bare 
Hamiltonian, Hb (with t = 0), is not quadratic in the 
boson operators. Therefore, one cannot use the standard 
diagrammatic techniques^ because the Wick's theorem 
does not hold. To make progress, we write the corrections 
to the Green's function using the cumulant expansion: 

({TMr)b\{Q)))=(TMr)b\m (12) 

+52 drij dT 2 Mjl(Tl-T2)Kij l (T,Ti,T2), 

jt Jo Jo 

where ((...)) denotes the Green's function which includes 
the perturbative corrections. In the Mott-insulating 



state, it is convenient to calculate the correlation function 
Kiji(T 1 T\,T2) in the second quantized representation: 

K ij i(T,n,T 2 ) =<T r 6 i (T)^(0)n i (Ti)n,(75)) (13) 
-(TAM^WKT^nMra)). 

Given that the on-site part of the boson Hubbard 
Hamiltonian conserves the number of bosons, the cor- 
relation functions above can be calculated exactly us- 
ing the particle-number eigenstates^i The terms in 
Kiji(r, Ti, T2) contributing to static and dynamic screen- 
ing are given by, 

Km (t, n , r 2 ) = 6 (t)0 (n ) e (r 2 ) 6 (r - n)6 (r - r 2 ) 
x [{5ij+8a)no(nQ+\)+5ii8u{na+\)]ex$ (-8E p t) 
+ e(-r)e(-T 1 )e(-T 2 )G(r 1 - r)e(r 2 - r) 
x [-(Sij + 8u)nl + 8ij8un ] exp (8E h r) + .... (14) 

It is important to note that Kiji(r, Tx, t 2 ) is irreducible 
and cannot be factored into the product of the bare 
Green's functions, as would have been possible if Wick's 
theorem were applicable. 

We now proceed to calculate the effects of the fermions 
by first approximating M g (f2 n ) in Eq. (O by the con- 

stant piece, M q (Cl n ) ~ -fg- (static approximation, see 
also Ref. [5.11]). By substituting the corresponding ex- 

u 2 

pression for Mji, M,;(ti-t 2 ) = -^-5ij&(T\ - r 2 ), into 
Eq. (d2), and carrying out the imaginary-time integrals, 
we find the following expression for the Green's function 
at zero frequency^, 



G((0) = - 



np + 1 
8E V 



2n )' 



2A8E„ 



n 
8E h 



1- 



'Fff 



(l-2noJ 



2A8E h 
(15) 



Alternatively, we could substitute the static, on-site form 
of Mij(r\ — r 2 ) directly into the action, Eq. ((5]), and cal- 
culate the Green's function exactly. It is easy to see that, 
in the static approximation, the mobile fermions simply 
renormalize fi and U of the bare boson Hubbard Hamil- 
tonian H B - U -> U - U FB /A and fj, -> fx + U FB /2A. 
The exact Green's function, thus, can simply be ob- 
tained by substituting these rcnormalizcd parameters in 
Eq. (fTTj) . After expanding the result to the second order 
in Ufb, the resulting expression exactly matches 21 that 
in Eq. (|T5]) . This validates the correctness of our pertur- 
bation theory. Using Eq. (JTTJJ) one can see that, in the 
static approximation, the fermions markedly shrink the 
area of the Mott-insulating lobes in the phase diagram 
(see Fig. [TJ. 

The static screening approximation for M q (f2„) does 
not, however, take into account the important retarda- 
tion effectsii and the spatially non-local nature of the 
interaction kernel in Eq. ([6]). By substituting the full 
expression for My(ri — r 2 ) into Eq. (fl2"j) . and doing the 
imaginary-time integrals as well as carrying out the sum- 
mation over j and Z, we obtain the following expression 
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for the boson Green's function at zero frequency, 



Gi(0) = 
n + l 



5E V 



1 + 



U 



FB 



ASE, 



R 



SE P 
AE F 



n 
SE h 



(16) 

U 2 fb nW 



Here we introduced the dimcnsionless function R(y)'. 



R(vY- 



kdk 

n JO 



(17) 



7+2/- ^2/ 2 +2/V / 2/ 2 -lsec 



The inset in Fig. [T] depicts the behavior of the monotonic 
function R{y) as a function of its argument. As follows 
from Eq. (|16[) . the importance of the fermion renormal- 
ization effects is determined by the ratio oi8E p /h and Ep. 
When the fermion density is small, i.e., SEp/^/Ep 1, 
the corrections to the Green's function are suppressed 
since R(y^> 1)— >0. In the opposite limit, SE p / h /Ep <C 1, 
the function R(y <§C 1)~1, and thus, for a given value of 
Ufb, the effects of the fermions on the bosons are more 
pronounced. Finally, using Eq. (Til)]) and Eq. (fTU)) , we cal- 
culate the phase diagram on the (fi — t) plane as shown in 
Fig. [1] We emphasize that the net effect of the fermions 
is to suppress the Mott-insulating lobes and enhance the 
superfluidity. 

The above result is consistent with numerical calcu- 
lation of Ref. |lfj| and is in disagreement with the con- 
clusions of Ref. [ljj. We note that the correctness of 
our formalism for the perturbative evaluation of the 
Green's function (in the static screening approximation) 
was confirmed independently, see the discussion after 
Eq. (fl~5)) . The generalization of the scheme to the dynam- 
ical screening is straightforward and amounts to only tak- 
ing the frequency and momentum integrals, mandated by 
Eq. (fl2)) . Thus, we are able to calculate the perturbative 



effects to the boson Hubbard model of an arbitrary time- 
and space-dependent interaction kernel. In contrast, it is 
not obvious that a spatially non-local interaction kernel, 
such as that in Eq. can be properly treated in the 
Weiss-like self-consistent mean-field theory employed in 
Ref. [ll|. Note also that the function M q (£l n ) is positive 
definite for all momenta and frequencies. Therefore, the 
net effect of the full interaction kernel is qualitatively sim- 
ilar to its constant piece (the static approximation), even 
though the latter significantly overestimates the suppres- 
sion of the insulating phase. Thus, our qualitative con- 
clusions should be valid for 3D systems as welljSi We note 
that the sign of the phase boundary shift can be predicted 
from the sign of the fermion density-density correlation 
function, while the magnitude of the corrections to the 
phase diagram depends on the microscopic details such 
as the ratio of 5E p / h and Ep as follows from Eqs. (fl6)) 
and (|17p . Finally, we emphasize that, near the degener- 
acy points, where the excitation energy SEp/^ is smaller 
than U FB /A, our perturbation theory breaks down, see 
Eqs. (fl~5|) and (fl6|) . Thus, the effect of fermions on the 
boson Hubbard phase diagram near these points is an 
open question. 

Conclusion. In summary, we develop a framework for 
carrying out the perturbation theory for the boson Hub- 
bard model, and use it to calculate the effects of a dilute 
gas of spin-polarized fermions weakly interacting with the 
bosons. The full theory captures both the static and the 
important dynamic effects of the fermions on the con- 
stituent bosons. We find that within single-band boson 
Hubbard model the net effect of the fermions is to in- 
herently suppress the Mott-insulating lobes and enhance 
the superfluid phase in the generic Bose-Hubbard phase 
diagram. 

We thank K. Yang, G. Refael, E. Dernier, T. Porto, 
T. Stanescu, E. Hwang, and C. W. Zhang for stimulating 
discussions. This work is supported by ARO-DARPA. 



1 M. Greiner, O. Mandel, T. Esslinger, Th. Hansen and 
I. Bloch, Nature 415 (2002). 

2 I. B. Spielman, W. D. Phillips, and J. V. Porto, Phys. Rev. 
Lett. 98, 080404 (2007). 

3 M. P. A. Fisher, P. B. Weichman, G. Grinstein, and D. S. 
Fisher, Phys. Rev. B 40, 546 (1989). 

4 S. Sachdev, Quantum Phase Transitions (Cambridge Uni- 
versity Press, Cambridge, 1999). 

5 H. P. Biichler and G. Blatter, Phys. Rev. Lett. 91, 130404 
(2003); Phys. Rev. A 69, 063603 (2004) 

6 M. Lewenstein, L. Santos, M. A. Baranov, and H. 
Fehrmann, Phys. Rev. Lett. 92, 050401 (2004). 

7 L. Mathey, D.-W. Wang, W. Hofstetter, M. D. Lukin, and 
E. Demler, Phys. Rev. Lett. 93, 120404 (2004); 

8 K. Sengupta, N. Dupuis and P. Majumdar, Phys. Rev. A 
75, 063625 (2007) 

9 A. Mering and M. Fleischhauer, Phys. Rev. A 77, 023601 
(2008) 

L. Pollet, C. Kollath, U. Schollwock, and M. Troyer, Phys. 



Rev. A 77, 023608 (2008). 

11 G. Refael and E. Demler, Phys. Rev. B 77, 144511 (2008). 

12 K. Giinter, T. Stoferle, H. Moritz, M. Kohl, and 
T. Esslinger, Phys. Rev. Lett. 96, 180402 (2006). 

13 S. Ospelkaus, C. Ospelkaus, O. Wille, M. Succo, P. Ernst, 
K. Sengstock, and K. Bongs, Phys. Rev. Lett. 96, 180403 
(2006). 

14 S. Chakravarty, G.-L. Ingold, S. Kivelson, and G. Zimanyi, 
Phys. Rev. B 37, 3283 (1988); S. Tewari, J. Toner, and S. 
Chakravarty, Phys. Rev. B 73, 064503 (2006). 

15 I. S. Beloborodov, A. V. Lopatin, V. M. Vinokur, and K. B. 
Efetov, Rev. Mod. Phys. 79, 469 (2007). 

16 S. Ospelkaus, C. Ospelkaus, L. Humbert, K. Sengstock, 
and K. Bongs, Phys. Rev. Lett. 97, 120403 (2006) 

17 M. Cramer, S. Ospelkaus, C. Ospelkaus, K. Bongs, K. Sen- 
gstock, and J. Eisert, Phys. Rev. Lett. 100, 140409 (2008). 

18 D.-S. Luhmann, K. Bongs, K. Sengstock, and 
D. Pfannkuche, Phys. Rev. Lett. 101, 050402 (2008) 

19 A. L. Fetter and J. D. Walecka, Quantum Theory of Many- 



Particle Systems (McGraw-Hill, New York, 1971). 

D. Dalidovich and M. P. Kennett. laxXiv:0711.1563l 

S. Tewari, R. Lutchyn and S. Das Sarma, (in preparation). 



Near the degeneracy points 3 where 8E v / h is of the 
der of Up B / A the perturbation theory breaks down. 



